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We consider the mixed inflaton-curvaton model in which quantum fluctuations of the curvaton 
field during inflation lead to a relatively large curvature perturbation spectrum at small scales. 
The predicted spectrum is characterized by the large positive tilt and the non-Gaussianity of the 
perturbations is also large. We obtained the constraints on the model parameters considering the 
process of primordial black hole (PBH) production in radiation era. The dependence of the produced 
' PBH's mass on model parameters is studied, as well as the form of predicted PBH mass spectra. 
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' Curvaton mechanism which has been suggested ^15 years ago flL's'l now is the object of intense study. It is assumed, 
in the standard implementation of the curvaton model, that not the inflaton field perturbations are responsible for the 
primordial density fluctuations and for the cosmic microwave background fluctuations, but instead the (isocurvature) 
perturbations of the curvaton field a. It is assumed that this curvaton field is subdominant during inflation but in 
post-inflationary epoch when Hubble constant becomes small, H m (where m is the curvaton mass), curvaton starts 
\ oscillating in its potential and behaves as nonrelativistic matter. The energy density of the curvaton decreases as 
Q-i, (a is the scale factor) whereas the energy density of radiation produced by the inflaton decay decreases as a~'^. 

Q As a result the curvaton energy density grows relative to radiation energy density until the curvaton contribution 
, becomes significant. If it happens before the curvaton decay one can say that curvaton mechanism is "effective" , 
C/3 ' in a sense that just the curvaton (rather than inflaton) field perturbations during inflation determine the resulting 
, ^, 1 (adiabatic) curvature perturbations at cosmological scales. 

In scenarios with the "effective" curvaton there is the strong constraint on a value of the curvaton mass: it must 
be much smaller than the Hubble constant during inflation. Hi, otherwise the primordial density perturbations have 
I> . too large spectral tilt. Moreover, if the ratio / Hf is not small, the coherent length of the curvaton field (i.e., 
0^ the characteristic size of the region inside of which the field is approximately homogenous) is also too small and, in 
, particular, smaller than the current horizon size. In the latter case, the primordial perturbation spectrum is strongly 
non-Gaussian, in contradiction with observations. 

The condition m^/Hf ^ 1 is too restrictive and prohibits an use, for a description of the curvaton, particle physics 
models predicting large ratios m?/Hf at inflation (e.g., some variants of supersymmetric theories). In this connection 
it is reasonable to consider also the mixed curvaton-inflaton scenarios 0, 0] in which the curvaton perturbations are 
additional to the usual perturbations produced by the inflaton. Combining two contributions, one can obtain the 
• • ' primordial perturbation spectrum which is in agreement with data at cosmological scales. At the same time, the 
^ : prediction for smaller scales may be quite unusual: the spectrum can be, e.g., very blue (i.e., the spectral tilt is large 
K> and positive) and, besides, the perturbations can be strongly non-Gaussian. In particular, large value of the tilt arises 
5_] , due to non-renormalizable and supergravity corrections to the Lagrangian of some supersymmetric theories inducing 
mass terms of the order \§Hl% . 

In most curvaton scenarios it is assumed that the curvaton field in the Universe is highly homogeneous and, as a 
result, the non-Gaussianity is small. Instead, we assume that, after the long inflationary expansion, the average value 
of the curvaton field is close to zero, and the local value of the field has a Gaussian probability distribution, variance 
of which is given by the formula [H, [13 

3Hf f H,V 1 
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Here, is the effective curvaton mass which differs from the true curvaton mass m The corresponding coherent 
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length is 

In Eqs. ([J) and ([2]), is the spectral tilt of the perturbation spectrum of the curvaton field, ta = dlnVa/dlnk. The 
assumption that (7 = will have real sense if the scale of interest, £ji = Ui/kji, will be larger than £c (both scales are 
calculated at the end of inflation). In our case, in is given by 

= T^e^. (3) 

Here, kend is the scale leaving the horizon at the end of inflation, is the scale factor at the end of inflation (and at 
the beginning of radiation era) , TV is a number of e-folds after the scale kji leaves the horizon. The condition in ^ £c 
leads to the inequality N ^ 1/ta- It means that if is not small {ta ~ 1), and the coherent length ic is small, one 
anticipates the blue curvature spectrum (the curvaton contribution) and large non-Gaussianity at small scales. In 
this case, the data at cosmological scales are described by the inflaton fluctuations only. 

Traditionally, predictions for the curvature perturbation spectrum in a region of small scales can be constrained 
with a help of primordial black holes (PBHs), because a large amplitude of the spectrum can lead to their production 



at radiation era. In the concrete case of the curvaton model, the idea was suggested in [16| and was considered, in 
more detail, in More recently, the PBH formation in an axion-like curvaton model was studied in [Tsj . 

In the present work we consider the PBH formation and PBH constraints in a special case when adiabatic perturba- 
tions at small scales are produced by curvaton only, resulting in a blue curvature spectrum and large non-Gaussianity. 
It means that in our case the typical size of the "curvaton domain" [l^ is relatively small, it is smaller than the 
horizon size at the moment of the formation of PBH with a given mass. 

The plan of the paper is as follows. In the next Section we derive the basic formula for the curvature perturbation 
spectrum used in the concrete calculations. In Sec. IIII[ the process of PBH production in our curvaton model is 
considered. The last Section contains the results of the calculation and conclusions. 



II. CURVATURE PERTURBATION SPECTRUM FORMULA 

The calculation of the curvature perturbation spectrum from the curvaton, V(^^ , is based on the expression for the 
nonlinear curvature perturbation on uniform density hypersurfaces (20l . [2l| , 

p(t,x) 

at,^) = 5N{t,^) + \ j (4) 
3 J p + p 

m 

where 6N is the perturbed expansion, p and p are the local density and the local pressure. From Eq. ^ we have 
formulas for the local curvaton and radiation densities on hypersurface on which curvaton decays (we use the sudden- 
decay approximation): 

C. = C+Jln^, (5) 

4 Pr 



C, = C+\\n^. (6) 

3 Pa 

Here, C is the total curvature perturbation after the decay. Besides, one has the expression for the local curvaton 
energy density on the hypersurface characterized by once the curvaton starts to oscillate 

p,=p,e'iC'-i^). (7) 

Using ([7]) , one can connect these perturbations with the field fluctuations during inflation. At the beginning of the 
oscillations one has 

Pae'^^-'^'-^ = \m^al,^, (8) 
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and the connection is given by the relation [2l| 



.9(0"*) = g{<7* + Sa^,). 



(9) 



Here we assume that the quantum fluctuations of the curvaton at the horizon exit during inflation are Gaussian, i.e., 
cr* = a* + (5 (7*, with no higher order (non- Gaussian) terms. 
We use the expansion [2l| 



where 



(10) 
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The final expression for Co- is (assuming that ^„ ^ 1) (see, e.g., [2: 
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where 
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The expression for i?^ is given by [1, [H, 
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(13) 



(14) 



where is the decay rate of the curvaton into radiation. It follows from (jl4p that, if <Tosc goes to zero, the linear 
term in Eq. (I12p vanishes. In this approximation, one has 



3 5^ 



(15) 



In what follows, we assume that the curvaton field is frozen during inflation, so we put cr* ~ ggi where CTg is the value 
of the field at the end of infiation. The connection of (Jcg-spectrum with the CTg-spectrum is |la] 



-,1/2 



1/2 



(16) 



and the spectrum of the curvaton field is [If 



(17) 



Here, Ninfi is the number of e-folds after the cosmological scales leave the horizon. The difference Ninfi — N is the 
number of "relevant" e-folds [16]. 

Finally, we obtain for the curvature spectrum the expression 



1.1/2^2 g'^ 1 Hf fk 



1/? 



(18) 



III. PBH PRODUCTION IN THE CURVATON MODEL 

We consider the model with the potential (see, e.g., the recent work [2^, where other examples of supergravity 
potentials are considered): 



(19) 
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The equation for the curvaton field a is 



a + 3H{t)a + V^=0. 



(20) 



The solution of this equation for the particular set of parameters is shown in Fig. [T] The calculation starts at moment 
t = corresponding to the end of inflation and the beginning of the radiation-dominated era (the reheating is assumed 
to be instant). During inflation, at t < 0, the solution of Ec^. (1^0]) is cr^ ~ exp{~taHit) [l^, and we use it to obtain 
the initial condition for ((^(7)) at t = 0. 

The spectral tilt ta is simply connected with a (if <C Hf): 



tn 



2ml 



3Hf 



(21) 



In this work we use, for concrete calculations, three values of a, a = 1, a = 0.6 and a — 0.38, so that the curvature 
spectrum is, in all these cases, blue, with the spectral index n = 1 + 2t„. 

The derivative g' is calculated numerically, from the solution of Eq. (PH)) . In our case, because the potential 
is quadratic, (Jose ^ fe, so g' — aosc/<Je- For the value of CTosciWe take cTosc = o'(iosc), with the moment of time when 
oscillations start, to 
is used) 



, which is determined by the condition [27[ (the first root of this Equation in the interval t > 



C^itosc) 



H{tosc)<j{tosc) 



= 1. 



(22) 



We are interested in the calculation of the PBH mass spectrum that are produced from the cosmological pertur- 
bations given by the spectrum of Eq. p^ . The general formula for the mass spectrum, in the Press-Schechter [26| 
formalism, is 12 8i 



UBHiMBH) = ( 



-1/3 



dP 
dR 



fhP. 



2/3^1/3 



(23) 



Here, pi,ai and Mi are initial (at time of reheating) energy density, scale factor and horizon mass, correspondingly, 
fh determines the relation between the produced PBH mass Mbh and horizon mass at this moment of time, Mh{t) 
(we assume in this work that //j w (1/3)"'^/^ — const): 



M. 



BH 



(24) 



and P{R) is the probability that in the region of comoving size R the smoothed value of curvature perturbation ^ 
will be larger than the PBH formation threshold value (^c- 



PiR) = / pciOdC- 



(25) 



Cc 



5 




In our case, when perturbations of C are quadratic in cr, one may write 

C = A(a2-(a2)), A>0 

so 

1 



PciO 



<7-nin ^ Cmin — A{a ), 



and 



Cmin — Cfnini^R^ 



\ J Vdk)W\kR)f 



1/2 



(26) 
(27) 

(28) 



where W{kR) is the Fourier transform of the window function, and we use a Gaussian one, W^{kR) ~ cxp(— fc^i?^), 
in this work. 



IV. RESULTS AND DISCUSSION 



The examples of the PBH mass spectra calculations are shown in Fig. [2] On the vertical axis, we show the 
combination of parameters giving, approximately, the value of Ppbh (energy density fraction of the Universe contained 
in PBHs at the time of their formation) [28j : 

1 /2 

Ppbh « ^ J nsHMl^dlnMBH « Ah^'^pi^AfJ^ ^bh{Mbh)\msh=m^^„- > (29) 

where M'^ff^s are the values of Mbh at maximums of curves in Fig. [2] For the curves shown in Fig. [21 the parameter 
Fo- was tuned so that the three curves corresponding to the different values of a gave approximately the same Ppbh- 

It is seen from Fig. [2] that for smaller values of a, the mass spectra become more wide. The low mass cut-off of 
the curves shown is determined by the fact that no PBHs are formed before the curvaton decays at t = tdec, so the 
minimal PBH mass is M™™ = fhMh{tdec)- 

For the constraining of the curvaton model parameters, we used the limits for (3pbh{Mbh) from the review work 
29|. Demanding that PBHs are not overproduced, i.e., the value of (3pbh{Mbh) does not exceed the available limits 
29j, one may obtain the corresponding constraints on the parameters of the considered cosmological model. Such 
constraints are shown in Fig. [3l 

One can see from the resulting Fig. [3]that, generally, PBH constraints are very weak. The forbidden region contains 
too small values of Trj/m (although the nucleosynthesis limit, ^ (lMeV)^/Mp, allows such values). The PBH 
constraint works only for very high values of Hubble constant during inflation. Hi > 10^° GeV, and for very large 
values of curvaton masses, m > (0.001 — 0.1)iJj. For other values of parameters, the spectrum amplitude, 7-"^^ is too 
small and cannot be constrained. One must note also that in the forbidden region the reheating temperatures are 
rather high {Tpn ^ vHiMp) and, in standard supersymmetric models, gravitinos are overproduced. 
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FIG. 3: The resulting constraints on the value of Vajm for the curvaton model considered in this paper. The shaded regions 
correspond to the sets of parameters that are prohibited by PBH overproduction. Upper curve corresponds to the case of 
m = 0.1//i, a = 0.38, middle curve - for m = 0.01-ffi, a = 0.38, lower one - for m = 0.01//i, a = \. Numbered labels show 
the PBH masses that are responsible for obtaining the constraints for the corresponding values oi Hi: 1 - Mbh ~ 10'^^ g; 2 - 
Mbh « 10^^ g; 3 - Mbh ~ 10^** g; 4 - Mbh « lO^^* g; 5 - Mbh ~ 10^^ g. 
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